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Some  Ray  Tracing  Problems  Related 
to  Circles  and  Ellipses 


1.  INTRODUCTION 

Practical  applications  of  optical  and  quasi-optical  methods  always  imply  solution  of  ray 
tracing  problems.  Such  problems  are  usually  very  easy  to  formulate,  but  analytical  solutions  are 
in  general  not  available.  They  are  therefore  normally  solved  numerically  and  In  practical 
applications  of  the  geometrical  optics  and  diffraction  theories  the  major  part  of  the  computer  time 
is  spent  solving  ray  tracing  problems. 

We  shall  consider  a  number  of  ray  tracing  problems  associated  with  circles  and  ellipses. 

Such  problems  arise  in  applications  involving  reflection  by  spheres  and  circular  and  elliptic 
cylinders  and  diffraction  by  circular  and  elliptic  edges.  We  shall  assume  that  these  bodies  are 
conductors  Imbedded  in  a  homogeneous  dielectric.  Thus  the  rays  we  are  to  consider  are  segments 
of  straight  lines  and  refracted  rays  are  excluded.  For  the  sake  of  convenience  we  shall  assume 
throughout  the  report  that  the  refractive  index  is  unity. 

Our  point  of  departure  will  be  a  classical  ray  tracing  problem,  namely  that  of  reflection  by  a 
spherical  mirror.  Although  not  well  known,  this  problem  was  originally  formulated  in  an 
ingenious  way  that  allowed  an  analytic  solution.  This  present  study  leads  to  an  alternative 
formulation  of  the  law  of  reflection,  which  directly  produces  the  simple  equation  for  the 
determination  of  the  reflection  point.  Furthermore,  an  alternative  ray  tracing  technique,  the  locus 
method,  is  introduced  and  applied  to  the  ray  tracing  problems  considered. 


(Received  for  publication  5  April  1989) 
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2.  ALHAZEN'S  MIRROR  PROBLEM 


2.1  The  Original  Formulation 

This  famous  problem  stems  from  the  Arabic  mathematician  Alhazen  (987-1038).  In  his 
Optics  the  problem  has  the  following  form:  Find  the  points  on  a  concave  spherical  mirror  at  which 
a  ray  of  light  coming  from  a  given  point  must  strike  in  order  to  be  reflected  to  another  given  point. 
The  beginnings  of  this  problem  are  found  in  Ptolemy's  works  on  optics,  and  after  Alhazen's 
masterly,  however  complicated,  discussion  of  it.  it  became  famous  in  Europe  on  account  of  the 
geometrical  difficulties  to  which  it  gave  rise  l*l. 

In  the  following  we  shall  review  Alhazen's  problem  in  the  form  given  by  Dorrie  W.In  Figure  1 . 
P^x,.  y,)  and  P2(x2.  y2)  are  the  two  given  points  while  R(x,  y)  is  a  reflection  point  on  the  circle  with 
radius  r.  For  later  use  we  also  introduce  polar  coordinates  for  these  points: 


ri  =Vx?+y? 

tan  a,  =  y!/x! 

r2  =VX2+Y2 

tan  a2  =y2/x2 

(1) 

r  =  ^x2+  y2 

tan  a  =  y/x. 

We  denote  the  angle  of  incidence  by  v,  and  the  angle  of  reflection  by  v2  and  furthermore 
introduce  the  angles  U!  and  u2  between  the  x-axis  and  the  rays  RPj  and  RP2,  respectively.  We 
then  have 

vi  =a-ut 
v2  =u2  -  a 

and  from  the  law  of  reflection  Vj  =  v2.  Introducing 

tan  a  =  — ;  tan  u,  =  — — — ;  tan  u,  =  - — —  io\ 

x  x-x,  2  x-x2  w 

into  the  equation 

tan  v,  =  tan  v,. 

(3) 

’Cajori.  F.  (1962)  A  History  of  Physics.  Dover  publications.  New  York,  pp.  21-23. 

2Dorrier.  H.  (1965)  A  Hundred  Great  Problems  of  Elementary  Mathematics.  Dover 
Publications.  New  York,  pp.  197-200. 
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Figure  1 .  Alhazen's  Mirror  Problem. 
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we  obtain  after  some  reductions 


xyi-yx,  xy2-yx2 

x2+y2-xx,-yy1  x2  +y2 -xx2 -yy2  ' 


(4) 


Now,  since  we  are  interested  only  in  the  reflection  points  on  the  circle,  we  may  substitute  r2 
for  x2  +  y2,  so  that  the  formula  can  be  written 

xyi-yxi  _  xy2-y*2 

r2-xxi-yyi  r2-xx2-yy2’  {5) 


Reduction  to  a  standard  form  yields 


Ax2  -  Ay2  +  2Cxy  +  r2(2Dx  +  2Ey)  =  0 


(6) 


where 


A  =  xjy2  +  x2y,  =  rtr2  sin(a,  +a2);  D  = -C.5(y,  +  y2) 

(7) 

C  =  yiy2 -XiX2  =-rir2c°s(a1 +a2);  E  =  0.5(X] +x2). 

Since  the  curves  given  by  Eq.  (6)  are  conic  sections  and  since  these  have  in  general  four 
points  of  intersection  with  a  circle,  it  may  be  concluded  that  there  are  at  most  four  reflection 
points  on  a  spherical  mirror. 

We  neglect  the  case  where  both  points  are  at  the  center  of  the  circle,  such  that  there  is  an 
infinity  of  reflection  points. 

This  result  was  obtained  by  Alhazen.  The  determination  of  the  reflection  points  involves 
solving  a  quartic  equation,  the  theory  of  which  was  not  known  to  Alhazen.  Thus,  as  we  shall  see, 
the  problem  had  to  wait  800  years  for  its  solution. 

2.2  Discussion  of  Alhazen' s  Curve 

The  conic  section  represented  by  Eq.  (6)  is  a  simple  hyperbola  with  asymptotes  that  are 
perpendicular  to  each  other.  Since  no  constant  term  is  present,  the  curve  contains  the  origin. 
From  the  eigenvalue  equation 


A-X  C 
C  -A-X 
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we  obtain  the  eigenvalues 


\  =  ±  Va2  +  C2  =  ±  r,r2 


(8) 


where  we  have  used  Eqs.  (7)  and  (1).  To  determine  the  center  (x,.,  yc)  of  the  hyperbola  we  solve  the 
linear  equations 

Axc  +  Cyc  +Dr2  =0 
Cxc  -  Ayc  +  Er2  =  0 


and  find,  again  using  Eq.  (7) 


(xc.yc)  = 


2 


(9) 


The  angle  to  by  which  we  rotate  the  coordinate  system  to  remove  the  xy-term  from  Eq.  (6)  is  found 
from 


which  yields  for  example 


tan  co  = 


C 

A  +  \ 


Inserting  X  =  r,r2  we  obtain,  after  some  reductions 


<o  =  0.5(a, 


(10) 


The  new  coordinate  system  C,  thus  has  its  origin  at  the  center  C  and  axes  which  arc  at 
45  degree  angles  with  the  angular  bisector  of  the  angle  defined  by  the  radii  through  P]  and  P2. 


It  may  be  shown  through  a  lengthy  calculation  that  the  equation  of  the  hyperbola  in  C  ]  is 


x2-y2=2k;  k  =  xcyc 


(ID 


where  (x^.,  yc)  are  the  coordinates  of  the  center  C  in  the  new  coordinate  system  C, .  given  by 


(Xc-yc)>“^r 4 


r2  +  rt  Oi+a2  r9-r,  .  a, +a 


-cos 


I  .  U]  -  u-2 

1  sin  — * - - 


(12) 


which  yields 


r4  r2-r2 

kr  ro  il  .  / 

=  7  =  JT2  Sin(a1+a2) 
o  rt  r2 


(13) 


Figure  2  shows  an  example  of  Alhazen's  hyperbola  for  r  =  10,  (xj,  y,)  =  (1,  6)  and  (x2,  yj  = 
(3.  -2).  The  axis  marked  x2  is  the  angular  bisector  of  the  angle  between  the  radii  through  Pj  and 
P2.  In  this  example  only  two  points  of  reflection  exist,  namely  P  and  Q,  only  one  branch  of  the 
hyperbola  intersects  the  circle. 
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REFLECTION  BY  CIRCLE:  INPUT  DATA 


x,  =  1.000  y,  =  6.000 

R,  =  6.083  a,  =  80.538 

x2  = 

R2  = 

3.000 

5.000 

Y2  = 

a2  = 

-4.000 

-53.130 

R  =  10.000 
v  =  13.704 

COEFFICIENTS  OF  POLYNOMIAL 
a=  14.000  c  =  -27.000 

d  = 

-1.000 

e  = 

2.000 

CONSTANTS  OF  HYPERBOLA 
xc  =  7.351  yc  =  0.108 

xc,  =  7.168  yc,  =  -1.637 

xm  = 

k  = 

2.000 

-11.730 

ym  = 

1.000 

ANALYTICAL  REFLECTION  POINTS  AT  FI  =  67.16  199.16 
NUMERICAL  REFLECTION  POINTS  AT  FI  =  67.16  199.16 

Figure  2.  Alhazen's  hyperbola  (dotted)  and  related  coordinate  systems. 
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2.3  An  Analytic  Solution  to  the  Mirror  Problem 

In  1861  the  Danish  mathematician  Tychsen  published  a  simple  analytic  solution  to  Alhazens 
problem  l3L  Probably  because  his  results  were  published  in  Danish,  they  never  became 
internationally  known.  His  solution  is  based  on  a  theorem  that  he  published  in  1851  Ml.  This 
theorem  states: 

When  the  real  coefficients  in  the  equation  of  the  fourth  degree 


x4  +ax3  +  bx  - 1  =  0 


(14) 


satisfy  the  condition 


a  +  bj3  ^a-bj3  |  t<Q  (15) 


the  equation  has  four  different  real  roots,  three  real  roots  or  two  real  roots,  respectively. 

In  order  to  utilize  this  theorem,  the  equation  of  the  hyperbola,  as  given  by  Eq.  (1 1)  is 
transformed  to  a  coordinate  system  C2  with  its  origin  at  the  center  of  the  circle  and  with  axes 
parallel  and  perpendicular  to  the  bisector  mentioned  earlier.  The  transformation  equation  is 


xl 
yJ2 


V2 

V2 

2 

2 

V2 

2 

2 

'cj. 


which  yields 


xy-ycx-xcy  =  0 


(17) 


in  the  coordinate  system  C2.  where  (x^.,  yc)  is  given  by  Eq.  (12).  Next,  to  derive  an  equation  of  the 
form  of  Eq.  (14).  a  new  variable  z  is  introduced,  and  we  rewrite  the  equation  of  the  circle  in  the 
form 


^Tychsen.  C.  (1859)  On  a  Mathematical  Investigation  of  a  Billiard  Problem.  (In  Danish). 
Matematisk  Tidsskrijt,  Copenhagen,  pp.  60-64. 

4Tychsen,  C.  (1851)  A  Remark  Concerning  an  Equation  of  the  Fourth  Degree.  (In  Danish). 
Matematisk  Tidsskrijt,  Copenhagen,  pp.  141-144. 
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(18) 


r  +  x  _  y 


y  r  -  x 


Solving  these  equations  for  x  and  y  we  find 


1 —  z*  2z 

x  =  r - 5-;  y  =  r - 5-. 

1  +  z2  J  1  +  z2 


(19) 


These  expressions,  when  inserted  into  Eq.  (17).  lead  to  the  fourth  degree  algebraic  equation 

4  „r  +  xr  3  „r-xf  ,  . 

2  "2— ^z3  +  2~— s-z-1  =  0.  ,20) 

y  c  y  c 


which  is  of  the  same  form  as  Eq.  (14).  Insertion  of  the  coefficients  in  the  condition  [Eq.  (15)1  then 
yields 


2  2  2 

x3+y3=r3  ^1) 


which  implies  that  the  problem  has  four  solutions  if  the  center  C  of  the  hyperbola  is  located 
inside  the  so-called  asteroid,  given  by 


2  2  2 
x3  +y3  =r3 


and  two  solutions,  if  C  is  outside.  This  is  Tychsen's  elegant  solution  to  Alhazen's  problem.  In 
Figure  3  the  asteroid  is  shown  for  the  same  values  of  the  parameters  as  were  used  in  Figure  2. 
Note  that  the  above  expression  for  the  asteroid  refers  to  the  coordinate  system  C2. 

It  is  worthwhile  to  note,  that  by  the  above  technique  the  problem  has  been  solved 
analytically,  since  the  solution  of  an  algebraic  quartic  equation  may  be  obtained  by  analytical 
methods.  This  is  due  to  the  fact  that  by  using  Alhazens  formulation,  no  square  roots  are 
introduced  in  Eq.  (4),  which  in  turn  implies  that  it  is  not  necessary  to  square  the  equation  in 
order  to  solve  it.  Squaring  immediately  raises  the  degree  of  the  equation  above  four  and  thus 
excludes  the  possibility  of  an  exact  solution.  In  the  next  section  we  shall  discuss  this  problem  in 
further  detail. 
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REFLECTION  BY  CIRCLE:  INPUT  DATA 


x,  =  1 .000  y,  =  6.000 

R,  =  6.083  a,  =  80.538 

x2  = 

r2  = 

3.000 

5.000 

y2  = 
a2  = 

-4.000 

-53.130 

R  =  10.000 
v  =  13.704 

COEFFICIENTS  OF  POLYNOMIAL 
a  =  14.000  c  =  -27.000 

d  = 

-1.000 

e  = 

2.000 

CONSTANTS  OF  HYPERBOLA 
xc  =  7.351  yc  =  0.108 

xc,  =  7.168  yc,  =  -1.637 

xm  = 

k  = 

2.000 

-11.730 

ym  = 

1.000 

ANALYTICAL  REFLECTION  POINTS  AT  FI  =  67.16  199.16 
NUMERICAL  REFLECTION  POINTS  AT  FI  =  67.16  199.16 


Figure  3.  Tychsen’s  Asteroid. 
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2.4  An  Alternative  Formulation  of  the  Law  of  Reflection 

In  most  textbooks  the  law  of  reflection  Is  usually  expressed  through  Snells  laws,  namely  that 
the  angle  of  incidence  and  the  angle  of  reflection  are  equal  and  that  the  planes  of  incidence  and 
reflection  coincide.  In  Figure  4,  S  is  a  source  point  and  F  a  field  point.  Using  the  reflection  point  R 
as  a  reference  the  position  vectors  for  S  and  F  may  be  written 


RS  =  S,  =  SjSj 
RF  =  S3  =  SiSi 


Furthermore  let 


N  =  Nn 


be  a  normal  vector  to  the  reflecting  surface  at  R.  The  law  of  reflection  then  is  usually  expressed 
by  the  formula 


n  x  s,  =s2xn 


(22) 


It  is  seen  that  square  roots  are  implicitly  involved  in  Eq.  (22)  and  in  order  to  solve  for  the 
coordinates  of  R  it  is  in  general  necessary  to  square  it. 

Another  formulation  of  the  law  of  reflection  may  be  based  directly  on  Fermat's  principle, 
which  states  that  the  optical  path  length  from  S  via  R  to  F  shall  be  stationary.  Referring  again  to 
Figure  4.  the  optical  path  length,  assuming  a  refractive  index  of  unity,  may  be  written 

s  =  s,  +s2 


where  s2  and  s2  are  given  by  square  roots.  These  are  preserved  after  the  differentiation  implied  by 
Fermat’s  principle  and  squaring  again  must  be  performed  to  solve  the  equation  for  the 
coordinates  of  R. 

It  seems  to  be  little  known  that  it  is  indeed  possible  to  formulate  the  law  of  reflection  without 
using  square  roots.  Referring  to  Figure  4.  we  rewrite  Eq.  (22)  in  the  form 


N  x  S[  _  S2  x  N 
NS!  "  S2N 


and  note  that 
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Figure  4.  The  Law  of  Reflection. 


N*S,  =  NS[  cos  v, 
N*S2  =  NS2  cos  v2 


where  Snell's  first  law  implies  that  vl  =  v2.  Inserting  In  the  above  equation  then  yields 


NxSt  _  S2  x  N 
N.S[  ~  S2.N 


which  is  the  expression  sought.  It  appears  that  it  contains  no  square  roots  and  thus  is 
particularly  useful  for  further  calculations. 

To  illustrate  the  use  of  the  formula,  we  apply  it  to  Alhazen's  problem.  We  have 


We  find 


NxS,  =(-xy!  +  yx,)z;  N*S!  =  x2  +  y2  - xxt  -  yy, 

NxS2  =(-xy2 +yx2)z:  N*S2  =  x2  +y2  -  xx2 -yy2 


Insertion  In  the  law  of  reflection  [Eq.  (23)1  then  immediately  yields  Alhazen's  equation  as 
given  by  Eq.  (4). 

If  the  problem  is  addressed  using  the  usual  law  of  reflection  one  obtains 


-xyi+yx,  _  -xy2+yx2 

•J(x,-x)2+(yi-y)2  t/(x2-x)2+(y2-y)2 

which  upon  squaring  yields 

(-xyi+yx,)2  =  (-xy2  +  yx2)2  (25) 

x2+y2+r12-2(xx,+yy1)  x2 +y2 +  r^ -2(xx2 +  yy2) 


This  is  an  equation  of  the  fourth  degree,  that  is.  one  degree  higher  than  Eq.  (4). 

Very  often  the  solution  is  based  directly  on  Fermat’s  principle.  Expressing  the  optical  path 
length  as  a  function  of  the  angle  a  in  Figure  1 ,  we  find 

s  =  s,+s2 


where 


s,  r2  +  r,2  -r(X]  cosa  +  y,  sin  a) 
s2  =  -^r2  +  r2  -r(x2  cosa  +  y2  sina) 


Differentiating  and  equating  the  derivative  to  zero  then  yields  the  equation 

-y^osa  +  xtsina  y2cosa-x2slna 

s,  s2 


(26) 
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which  is  Identical  to  Eq.  (24).  The  inexpediency  of  Eq.  (26).  in  this  context,  appears  when  a 
solution  for  cos  a  or  sin  a  is  attempted.  This  leads  to  an  equation  of  the  sixth  degree  and  as  it 
well  known  such  an  equation  cannot  in  general  be  solved  analytically. 


2.5  The  Locus  Solution  to  the  Mirror  Problem 

In  this  section  we  shall  investigate  an  alternative  solution  to  Alhazen’s  problem. 
Resubstituting  x2  +  y2  for  r2  in  Eq.  (5)  it  is  turned  into  an  equation  for  the  locus  of  the  points  for 
which  Vj  =  v2,  (see  Figure  1).  The  equation  becomes 


2Dx3  +  2Ex2y  +  2Dxy2  +  2Ey3  +  Ax2  +  2Cxy  -  Ay  2  =  0  127) 

that  is.  an  algebraic  equation  of  the  third  degree.  It  appears  that  the  origin  is  contained  in  the 
curve  and  that  a  maximum  of  three  y-values  corresponds  to  each  value  of  x.  The  equation  may 
easily  be  parametrized  by  introducing 

y  =  xt ;  t  =  tan  a 


which  yields 


x  =  y0 


t2  +  2at  - 1 
(t-t0)(l  +  t2) 


(29) 


where 


v  =  A  Xiy2+x2y, 

0  2E  X[  +  x2 

a  C.  xix2-y,y2 
A  x,y2  +  x2y, 

t  -  p_yi+y2 
0  E  xt  +  x2 


This  surprisingly  simple  equation  shows  that  the  curve  has  a  double  point  at  the  origin.  The 
corresponding  values  of  the  parameter  t  are  found  from  the  equation 


t2  +2at  -1  =  0 


(30) 
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with  the  solutions 


t,  =  tan 


a.  +a9 


and 


to  = -cot 


a,  +a2 


Since  at  the  origin  t  equals  the  slope  of  the  curve  it  is  seen  that  the  curve  intersects  itself  at 
right  angles.  Furthermore  it  is  immediately  found  that: 

(1)  The  curve  intersects  the  x-axis  at 


x  =  —  =  —1-^?- for  t  =  0 
to  Yl+Y2 


(2)  The  curve  intersects  the  y-axis  at 


y*yo 


xiy2  +  x2y, 

x,  +x2 


for  t  = 


(3)  The  line 


y  =  t0x  = 


yi+y2 
x,  +  x2 


x 


is  an  asymptote.  The  asymptote  is  the  line  through  the  origin  and  the  midpoint  M 
between  P,  and  P2. 

(4)  The  points  Pj  and  P2  both  lie  on  the  curve.  This  takes  some  calculations. 

(5)  Since  the  curve  has  only  one  asymptote  and  only  one  double  point,  it  must  have  a 
closed  loop. 

To  determine  the  value  of  t  corresponding  to  the  peak  of  the  loop  we  note  that 
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r  =  |x|Vl  +  t2 


t2  +  2at  -  1 

y°  (t-to)VTTt2 


(31) 


Equating  the  derivative  to  zero  we  obtain  the  third  degree  equation 


(2a  + 10  )t3  -  3t2  +  3t0t  -1  =  0 


(32) 


for  t  =  tjnax.  This  equation  can  be  solved  analytically  and  may  be  determined  using  Eq.  (31). 
We  have  then  found  the  following  simple  alternative  rule:  The  mirror  problem  has  two  solutions  if 
r  >  rmax  and  f°ur  solutions  if  r  <  rmax. 

In  Figure  5  the  locus  is  shown  for  (x,.  yj  =  (5,  1)  and  (x2.  yj  =  (-4.5,  1).  The  asymptote  is 
given  by  y  =  4x  and  the  intersection  points  with  the  axes  are  (-0.25.  0)  and  (0.  -1).  It  appears  that 
since  the  radius  of  the  circle  is  smaller  than  the  height  of  the  loop  there  are  four  reflection  points. 
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REFLECTION  BY  CIRCLE:  INPUT  DAIm 
x,  =  5.000  y,  =  1.000 

R,=  5.099  a,  =  11.310 

33  * 

II  " 

-4.500 

4.610 

y2  = 

a2  = 

1.000 

167.471 

R  = 

v  = 

10.000 

89.391 

COEFFICIENTS  OF  POLYNOMIAL 
a  =  0.500  c  =  23.500 

d  = 

-1.000 

e  = 

0.250 

CONSTANTS  OF  HYPERBOLA 
xc  =  -0.973  yc  =  4.276 

xc,  =  4.265  yc,  =  1.018 

xm  = 

k  = 

0.250 

4.343 

ym  = 

1.000 

ANALYTICAL  REFLECTION  POINTS  AT  FI 
NUMERICAL  REFLECTION  POINTS  AT  FI 

=  265.30 
=  21.98 

150.53 

99.75 

99.75 

150.53 

21.98 

265.30 

Figure  5.  The  locus  (full)  and  Alhazen's  Hyperbola  (dotted)  for  the  Mirror  Problem. 
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3.  RAT  TRACING  OF  REFLECTED  RAYS 


3.1  The  Sphere  Problem 

The  point  of  reflection  on  the  convex  side  of  a  sphere  may  be  determined  analytically  by  a 
slight  modification  of  Alhazen's  and  Tychsen's  theory.  Since  the  reflection  point  must  lie  in  the 
plane  determined  by  the  center  of  the  sphere,  the  source  point,  and  the  Reid  point,  the  problem  is 
immediately  reduced  to  a  plane  problem.  Thus,  as  shown  in  Figure  6,  the  points  Pi  and  P2are 
now  outside  the  circle,  but  by  extending  the  line  segments  P)R  and  P2R  to  intersect  the  x-axis  we 
may  again  introduce  the  angles  U!  and  u2  and  use  Alhazen's  technique  to  determine  the 
hyperbola  given  by  Eq.  (6).  Thus,  formally,  the  problem  again  has  in  general  four  solutions.  As 
illustrated  in  Figure  6,  one  of  these  corresponds  to  the  actual  point  of  reflection  R]  and  another, 
R2,  to  the  reflection  by  the  concave  side  of  the  circle.  The  two  last  solutions,  R3  and  R^  are  false 
and  correspond  to  non-physical  solutions  where  tan  v,  =  tan  v2. 

To  determine  the  actual  point  of  reflection.  Eq.  (20)  is  solved,  again  by  standard  analytical 
techniques.  The  angular  coordinate  a  of  the  reflection  point  then  is  given  by 


1  - z2  2z 

cosa  = - tt',  sma  = - y 

1  +  z2  1  +  z2 


(33) 


The  proper  root  may  be  determined  by  the  requirement  that  it  should  be  in  the  smallest 
angular  region  determined  by  the  circle  and  the  two  points  Pj  and  P2. 


3.2  The  Cylinder  Problem 

3.2.1  THE  BASIC  EQUATION 

In  this  section  the  theory  outlined  in  Section  2  is  applied  to  the  ellipse.  The  ellipse  ray 
tracing  problem  arises  in  connection  with  the  determination  of  the  point  of  reflection  on  either 
side  of  a  circular  or  an  elliptic  cylindrical  reflector. 

To  establish  Alhazen's  equation  for  the  ellipse  we  apply  the  reflection  law  as  given  by  Eq. 
(23).  Assuming  the  ellipse  is  given  by 

b2x2  +a2y2  =  a2b2. 


where  a  and  b  are  the  semi  axes,  a  normal  may  be  written 


From  Eq.  (23)  we  then  find  the  expression 
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REFLECTION  BY  CIRCLE:  INPUT  DATA 


x,  =  8.000  y,  =  12.000 

R,  =  14.422  a,  =  56.310 

x2  = 

r2  = 

13.000 

13.342 

y2  = 

a2  = 

3.000 

12.995 

R  = 

v  = 

10.000 

34.652 

COEFFICIENTS  OF  POLYNOMIAL 
a  =  180.000  c  =  -68.000 

d  = 

-7.500 

e  = 

10.500 

CONSTANTS  OF  HYPERBOLA 
xc  =  5.575  yc  =  3.727 

xc,  =  6.705  yc,  =  -0.104 

xm  = 

k  = 

10.500 

-0.695 

ym  = 

7.500 

ANALYTICAL  REFLECTION  POINTS  AT  FI 

=  83.66 

32.35 

347.51 

215.01 

NUMERICAL  REFLECTION  POINTS  AT  FI  =  32.85  21 5.01 

Figure  6.  Reflection  by  sphere.  Alhazen's  curve  is  shown  dotted,  and  the  Locus  by  a  solid  line. 
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(35) 


(b2  -a2)xy-b2y,x  +  a2x1y  (b2  -a2)xy  -  b2y2x  +  a2x2y 
b2x2  +  a2y2  -  b2x,x  -  a2y,y  +  b2x2  +  a2y2  -  b2x2x  -  a2y2y 


To  obtain  Alhazen's  equation  corresponding  to  the  ellipse,  we  introduce  the  parameter 

c4  =  a2b2 

which  we  substitute  for  b2x2  =  a2y2  in  the  above  equation.  This  yields 

(b2 -a2)xy-b2y1x  +  a2x1y  (b2 -a2)xy-b2y2x  +  a2x2y  ^  (36) 

4  i  9  2  4.9  9  tt, 

c  -b  x,x-a  yj  c  -b  x2x-a  y2y 

which  may  be  written 

a!X2y  +  a2xy2  +  a3x2  +  a4xy  +  a5y2  +  c4(a6xy +  a7x  +  a8y)  =  0  (37) 


where 


a,  =2b2(a2-b2)E  a2  =  2a2(b2 -a2)D 

a3  =  b4A  a4=2a2b2C 

a5  =  -a4A  a6  =  2(b2  -  a2)  (38) 

a7  =  2b2D  a8  =  2a2E. 


It  appears  that  due  to  the  more  complicated  geometry  of  the  ellipse  the  degree  of  Alhazen's 
equation  has  been  raised  by  one,  that  is,  to  the  third  degree,  although  the  terms  containing  x3 
and  y3  are  not  present.  Furthermore,  an  additional  term  agxy  is  introduced,  but  for  a  =  b  =  r  the 
equation  reduces  to  Eq.  (6)  as  should  be  expected. 


3.2.2  THE  ALHAZEN  CURVE 

To  obtain  an  expression  for  Alhazen's  curve  we  rearrange  Eq.  (37)  to  yield 
b^y  +  b2xy2  +  b3x2  +  b4xy  +  b5y2  +  b6x  +  b7y  =  0 
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where 


b,  =b2(a2-b2)(x, +x2)  b2  =  a2(a2  -  b2)(y,  +y2) 

b3  =  b4(x,y2  +  x2y;)  b4  =2a2b2(y!y2  -x,x2  +b2  -a2) 

b5=-a4(x,y2  +  x2y,)  b6  = -a2b4(y! +y2)  (40) 

b7  =  a4b2(x,  +x2). 

and  from  which  it  appears  that  the  curve  contains  the  origin.  Since  Eq.  (39)  is  a  second  degree 
equation  in  y: 

(b2x  +  b5  )y2  +  (b,x2  +  b4x  +  b7  )y  +  (b3x2  +  b6x)  =  0.  (41) 


in  general  two  values  of  y  correspond  to  each  value  of  x,  and  two  intervals  of  x  may  exist  where 
Alhazen's  curve  is  not  defined.  Furthermore,  the  curve  intersects  the  x-axis  at  (xq.  0).  where 


__be  -2  yi+y2 


xn  =-~  =  a 


b3  x2y,+xiy2 


and  it  appears  that  in  general  there  is  a  vertical  asymptote  at 


b5  _  a2(x|y2  +  x2yt) 


b2  (a2  -b2)(y,  +y2) 


Similarly  the  equation  for  x 


(b,y  +  b3)x2  +(b2y2  +  b4y  +  b6)x  +(b5y2  +  b7y)  =  0, 


shows  that  the  curve  intersects  the  y-axis  at  (0,  yj,  where 


b7  .2  x,  +  x, 
y0=-_2.  =  b2 - 1 - 2— 

b5  x2y,+x,y2 
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and  that  in  general  a  horizontal  asymptote  exists  lor 


bg  =  b2(x1y  +  x2y1)  (46) 

bi  (a2 -b2)(x,  +  x2) 


Parametrizing  the  equation  by  y  =  xt  we  obtain 

x[(b2t2  +  b,t)x2  +  (b5t2  +  b4t  +  b3  )x  +  (b7t  +  b6)]  =  0  (47) 

from  which  it  appears  that  in  a  given  direction  there  are  in  general  two  curve  points  although 
sectors  may  exist  where  Alhazen’s  curve  is  not  defined.  Furthermore,  the  equation  shows  that  an 
oblique  asymptote  appears  in  the  direction  given  by 


b,  b2(x,  +  x2) 
b2  a2(y,+y2)' 


(48) 


This  asymptote  passes  through  the  point  (0,  q).  where 


Thus  it  is  possible  to  gain  a  broad  idea  of  the  complicated  nature  of  Alhazen's  curve,  but  an 
analytical  solution  to  the  problem  seems  out  of  reach  in  the  general  case. 

Figure  7  shows  an  example  of  Alhazen's  curve  for  the  ellipse,  namely  for  (x,,  yj  =  (5.0,  0.8), 
(x2.  y2)  =  (-4.5,  0.8).  a  =  10  and  b  =  8.  Although  these  values  are  similar  to  those  used  for  the 
circle  in  Figure  5.  the  general  picture  is  entirely  different.  It  even  appears  that  six  reflection 
points  may  exist  on  an  ellipse.  In  the  next  section  we  shall  look  further  into  this  problem. 


3.2.3  THE  LOCUS 

We  next  consider  the  locus  given  by  Eq.  (37)  with 
c4  =  b2x2 +  a2y2. 

Inserting  this  and  rearranging  the  equation  we  obtain 
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Asymptotes: 
x  =  0.69 
y  =  -1.42 
y  =  -0.2x  +  29.9 


-100 


REFLECTION  BY  ELLIPSE:  INPUT  DATA 


x,  =  5.000 

y,  =  0.800 

x2=  -4.500 

y2  = 

0.800 

a  =  10.000 

b  =  8.000 

COEFFICIENTS  OF  POLYNOMIAL 

a  =  0.40000E+00 

c  =  0.23140E+02 

d  =  -0.80000E+00 

e  = 

0.25000E+00 

CONSTANTS  OF  PSEUDO-LOCUS  CURVE 

d,  =  -0.65536E+04 

d2  =  0.43520E+04 

d3=  -0.44800E+04 

d4  = 

0.50000E+04 

d5  =  0.16384E+04 

d6=  -0.1 6461 E+06 

d7=  -0.40000E+04 

t,  =  -0.41162E+02 

t2  =  0.9951 2E-02 

to=  0.11405E+01 

x0  =  0.25000E+00 

y0  =  0.80000E+00 

Q  =  -0.1 4021 E+02 

NUMERICAL  REFLECTION  POINTS  AT  FI  = 

84.68  200.16 

279.14 

331.75 

Figure  7.  Reflections  by  an  ellipse.  Alhazen's  curve. 
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c1x3y  +  c2xy3  +  c3x3  +c4x2y  +  c5xy2  +  c6y3  +c7x2  +  c8xy  +  c9y2  =0. 


(49) 


where 


c,  = -2b2(a2  -  b2)  c2  = -2a2(a2  -  b2] 

c3  =  -b4 (y,  +  y2 )  c4  =  b2(2a2  -  b2 )(x,  +  x2 ) 

c5  =a2(a2  -2b2)(y,  +y2)  c6  =a4(x,  +x2) 

c7  =  b4(x,y2+x2y1)  c8  =  2a2b2(y,y2  - x,x2) 

c9  =  -a4(x,y2  +  x2y1). 


(50) 


Again  due  to  the  ellipse  geometry  the  degree  of  the  locus  has  been  raised  by  one,  that  is,  to 
the  fourth  degree  although  terms  containing  x4  and  y4  are  not  present.  The  origin  still  is  on  the 
curve,  but  in  the  present  case  this  point  is  a  double  point.  Since  Eq.  (49)  is  a  third  order  equation 
in  y: 


(c2x  +  Cg  )y3  +  (c5x  +  c9  )y2  +  (c,x3  +  c4x2  +  c8x)y  +  (c3x3  +  c7x2 )  =  0. 


(51) 


a  minimum  of  one  and  a  maximum  of  three  values  of  y  correspond  to  each  value  of  x.  Apart  from 
the  origin,  the  curve  intersects  the  x-axis  at  (Xq,  0).  where 


.  c?  _  Xiy2+x2yt 
c3  Y1+Y2 


(52) 


and  it  appears  that  in  general  there  is  a  vertical  asymptote  at 


£jj_  _  a2  (X1  +X2) 

C2  2(a2  -  b2) 


Similarly  the  equation  for  x, 


(c,y  +  c3  )x3  +  (c4y  +  c7  )x2  +  (c2y3  +  c5y2  +  c8y)x  +  (c6y3  +  c9y2 )  =  0. 


(53) 


24 


indicates  an  intersection  of  the  y-axls  at  (0,  yo).  where 


c9  x,y2  +  x2yt 
c6  x,+x2 


and  a  horizontal  asymptote  at 


c.i  t>2(yi+y2) 
ci  2(a2  -  b2) 


Parametrizing  the  equation  by  y  =  xt  we  obtain 


x2[(c2t3  +  cit)x2+(c6t3+c5t2+c4t  +  c3)  +  (c9t2  +  c8t  +  c7)j  =  0. 


Thus  in  general  there  are  two  curve  points  in  each  direction  t.  Since  the  equation 


(54) 


(55) 


(56) 


t2^ 


b2 


has  no  real  solution,  no  oblique  asymptote  is  present.  Thus  in  the  ellipse  case  the  locus  is 
somewhat  simpler  than  Alhazen's  curve.  In  Figure  8  the  locus  is  shown  for  the  same  points  P, 
and  P2  as  in  Figure  6.  It  appears  that  the  locus  has  a  loop  as  in  the  circle  case,  but  also  that  four 
reflection  points  exist  although  the  ellipse  does  not  intersect  the  loop  but  rather  the  additional 
branch  of  the  curve  now  present. 

We  next  return  to  Eq.  (49)  and  note  that  the  coefficients  c,  and  c2  do  not  depend  on  the 
coordinates  of  P,  and  P2 .  The  two  first  terms  of  Eq.  (49)  thus  may  be  written 


c,x3y +  c2xy3  =  -2(a2  -  b2)(b2x2  +a2y2)xy 


(57) 


We  may  now  use  Alhazen’s  argument  that  we  are  interested  only  in  the  points  of  intersection 
between  the  locus  and  the  ellipse  and  introduce  Eq.  (34)  for  the  ellipse 

i-2  2  2  2  2 1-2  4 

b  x  +a  y  =a  b  =c 
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REFLECTION  BY  ELLIPSE:  INPUT  DATA 

x,  =  5.000  y,  =  0.800  x2  =  -4.500  y2  =  0.800 

a  =  10.000  b  =  8.000 


COEFFICIENTS  OF  POLYNOMIAL 
a=  0.40000E+00  c  =  0.231 40E+02 

d  =  -0.80000E+00 

e  =  0.25000E+00 

CONSTANTS  OF  PSEUDO-LOCUS  CURVE 
d,  =  -0.65536E+04  d2  =  0.43520E+04 

d3=  -0.44800E+04 

d4=  0.50000E+04 

d5  =  0.16384E+04 

d6=  -0.1 6461 E+06 

d7=  -0.40000E+04 

t,  =  -0.41162E+02 

t2  =  0.9951 2E-02 

t0=  0.11405E+01 

x0  =  0.25000E+00 

y0=  0.80000E+00 

Q  =  -0.1 4021 E+02 

NUMERICAL  REFLECTION  POINTS  AT  FI  =  84.68  200.16  279.14  331.75 


Figure  8.  The  locus  in  the  ellipse  case. 


26 


and  obtain 


c,x3y  +  c2xy3  =  -2(a2  -  b2)c4xy. 


Introducing  this  in  Eq.  (49)  yields  what  we  shall  call  the  pseudo-locus.  Rearranging  the 
terms  we  find 


d1x3  +  d2x2y  +  d3xy2  +d4y3  +d5x2  +  d6xy  +  d7y2  =0 


where 


dj  =  -b4(y,  +y2)  d2  =  b2(2a2  -  b2)(x,  +  x2) 

d3  =  a2(a2-2b2)(y,+y2)  d4  =a4(x,+x2) 

d5  =  b4(x,y2  +x2y,):  d6  =2a2b2(y,y2  -  XjX2  -  a2  +  b2) 

d7  =  -a4(x,y2+x2y1). 


A  comparison  with  Eq.  (27)  for  the  locus  in  the  circle  case  shows  that  the  two  equations  are 
formally  identical.  In  the  circle  case,  however,  d,  =  d3,  d2  =  d4  and  d5  =  -d7.  so  the  geometry  of  the 
ellipse  disturbs  the  symmetry  of  the  locus  equation  and  prevents  a  thorough  analysis  of  the  curve 
as  in  the  circle  case. 

Rearranging  the  terms  in  descending  powers  of  y  yields 

d4y3  +  (d3x  +  d7  )y2  +  (d2x2  +  d6x)y  +  d,  x3  +  d5x2  =  0.  (6  n 


which  shows  that  in  general  there  are  three  values  of  y  for  each  value  of  x.  The  curve  intersects 
the  x-axis  at  the  origin  and  at  the  point  (xq,  0).  where 


di  y!+y2 


and  it  has  no  vertical  asymptotes.  The  analogous  equation  in  x. 


d,x3  +(d2y  +  d5)x2  +(d3y2  +  d6y)x  +  d4y3  +  d7y2  =  0. 
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shows  that  in  general  three  values  of  x  correspond  to  each  value  of  y.  Furthermore  the  curve 
intersects  the  y-axis  at  (0,  yo),  where 


d7  xiy2  +  x2yt 

0  d4  x,  +  x2 


(64) 


and  it  has  no  horizontal  asymptotes.  Parametrizing  the  equation  by  y  =  xt  and  solving  for  x,  we  get 


d7t2  +d6t  +  d5 
d4t3  +d3t2  +  d2t  +  d. 


(65) 


This  equation  reveals  several  interesting  facts.  First,  it  shows  that  there  is  only  one 
curvepoint  in  any  direction.  Next,  since  the  discriminator  of  the  numerator  is  always  positive,  the 
curve  always  passes  through  the  origin  and  this  point  is  a  double  point.  Furthermore,  at  the 
origin  the  tangents  are  no  longer  at  right  angles  to  each  other.  Finally,  it  appears  that  since  the 
equation 


d4t3  +  d3t2  +d2t  +  d,  =0 


always  has  one  real  solution  t  =  ^  there  will  always  be  at  least  one  oblique  asymptote  with  a  slope 
determined  by  t0.  It  is  easily  shown  that  this  asymptote  passes  through  the  point  (0,  q).  where 

q  -  ^7^0  +  dfitp  +  ^5 

3d4t2  +2d3t0  +d2 


Since  the  presence  of  more  than  one  oblique  asymptote  is  in  general  incompatible  with  the 
fact  that  the  curve  has  only  one  point  in  each  direction  we  conclude  that  the  pseudo-locus  has  a 
closed  loop.  This  implies  that  there  is  a  maximum  of  four  intersection  points  wtth  the  ellipse. 
Furthermore  we  may  conclude  that  the  two  branches  of  Alhazens  curve  not  passing  through  the 
origin  cannot  both  intersect  the  ellipse  (see  Figure  7).  Figure  9  shows  the  pseudo-locus  for  the 
same  points  P!  and  P2  as  in  Figures  7  and  8. 

The  above  results  suggest  an  attempt  to  solve  the  problem  by  using  the  multiplicative 
mapping  that  maps  the  ellipse  onto  the  circle  x2  +  y2  =  a2. 

Application  of  this  mapping  to  the  pseudo-locus  results  in  a  geometrical  configuration  very 
similar  to  that  of  Figure  5  for  the  circle.  It  turns  out.  however,  that  the  equation  analogous  to  Eq. 
(27)  is  an  algebraic  equation  of  a  degree  higher  than  four.  It  is  also  worthwhile  to  note  that 
although  the  points  of  intersection  between  the  ellipse  and  the  pseudo-locus  represent  reflection 
points  on  the  ellipse,  the  points  of  intersection  between  the  mapped  pseudo-locus  and  the  circle 
x2  +  y2  =  a2  do  not  in  general  correspond  to  reflection  points  on  the  circle.  Neither  the  original 
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REFLECTION  BY  ELLIPSE:  INPUT  DATA 

x,  »  5.000  y,  =  0.800  x2  =  -4.500  y2  =  0.800 

a  =  10.000  b  =  8.000 

COEFFICIENTS  OF  POLYNOMIAL 

a  =0.40000E+00  c=  0.23140E+02  d  = -0.80000E+00  e=  0.25000E+00 
CONSTANTS  OF  PSEUDO-LOCUS  CURVE 

d,  .  -0.65536E+04  d2  =  0.43520E+04  d3  =  -0.44800E+04  d4  =  0.50000E+04 

d5  »  0.16384E+04  d6=  -0.1 6461 E+06  d7=  -0.40000E+04 

t,=  -0.41162E+02  t2=  0.9951 2E-02  to  =  0.11405E+01 

x0=  0.25000E+00  y0  *  0.80000E+00  Q=  -0.14021E+02 

NUMERICAL  REFLECTION  POINTS  AT  FI  =  84.68  200.16  279.14  331.75 


Figure  9.  Reflection  by  an  Ellipse.  Alhazen's  curve  (dotted),  and  the  Pseudo-locus  (solid). 
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points  P[  and  P2  nor  the  points  obtained  by  applying  the  mapping  to  these  points  are  reflection 
points.  On  the  other  hand  the  mapping  may  be  used  for  an  a  priori  determination  of  the  number 
of  reflection  points. 


4.  TRACING  OF  DIFFRACTED  RATS 
4.1  The  Basic  Equation 

The  generalization  of  Alhazen’s  problem  to  the  ray  tracing  problem  associated  with 
diffraction  by  a  circular  edge  is  based  on  the  similarities  between  the  laws  of  reflecUon  and 
diffraction.  These  laws  may  be  written 


s1xn  +  s2xn  =  0  and  S[*t  +  s2*t  =  0 


where  s,  is  a  unit  vector  pointing  from  the  reflection-  or  the  diffraction  point-  to  P,.  S2  points 
from  the  same  point  to  P2,  n  is  a  unit  vector  normal  to  the  circle  (in  the  plane  of  the  circle),  while 

A 

t  is  a  unit  vector  tangent  to  the  circle.  The  reflection  law  states  that  the  angles  of  incidence  and 
reflection.  v(  and  v2,  are  equal  and  that  the  planes  of  incidence  and  reflection  coincide.  The 
diffraction  law  states  that  the  angles  of  incidence  and  diffraction,  v,  and  v2  are  equal  and  that  the 
planes  of  incidence  and  diffraction  do  not  in  general  coincide.  Since  Alhazen's  theory  relies  on  the 
angular  relation  v,  =  v2  and  not  explicity  on  the  vectorial  nature  of  the  reflection  law.  it  may  be 
generalized  to  include  the  diffraction  problem  also. 

In  Figure  10.  the  circular  edge  is  given  by  x2  +  y2  =  r2  and  the  point  D,  with  coordinates  (xD, 
Yd-  zd)  =  (r  cos  a,  r  sin  a.  0),  is  assumed  to  be  a  point  of  diffraction.  At  this  point  we  introduce  a 
rectangular  coordinate  system  (x,,  y,.  z^  with  die  z,-axis  tangential  to  the  circle  and  the  yraxis 
perpendicular  to  the  xy-plane.  The  zraxis  thus  is  the  axis  of  the  Keller  cone  with  its  apex  at  D. 
The  cone  intersects  the  xy-plane  in  the  line  1  (and  1 ,)  and  we  denote  the  angle  between  this  line 
and  the  line  OD  by  v.  The  apex  angle  of  the  cone  then  is  2(90  -  v)  and  the  equation  of  the  cone 
becomes 


x?  +  Y?  -cot2vz,  =0 


(66) 


To  find  the  equation  of  the  Keller  cone  in  the  xyz-system  we  apply  the  transformation 


*1 

cos  a 

sin  a 

0 

x-xD 

Yi 

■  = 

0 

0 

1 

■  +  ■ 

y-yD 

zi 

-sin  a 

cos  a 

0 

z 

(67) 


and  find  after  some  reductions 
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z 


Figure  10.  The  Keller  Cone. 


tan2v  =  ■ 


(yxD-xyD) 


(r2  -  xxD  -  yyD)  +  r2z2 


(68) 


We  next  Introduce  the  two  points  Pi(xlt  yt,  z,)  and  P2(x2,  y2.  Z2)  and  the  corresponding  angles 
v,  and  v2.  Since  these  points  must  lie  on  the  cone  [Eq.  (66))  we  find  from  the  law  of  diffraction 
that  v,  =  v2,  or 


tan2v,  =  tan2v2 


(69) 


and  that 
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(y,xD  -x,yD) 


(y2xD-x2yD) 


(r2  -x,xD  -y^u)  +  r2z?  (r2  -x2xD  -y2yD)  +  r2z 


(70) 


It  is  noted,  that  for  zx  =  ^  =  0  the  formula  is  reduced  to  Eq.  (4),  however  squared. 

To  investigate  the  formula  derived  in  more  detail,  we  introduce  the  following  functions 

fi(x-y)  =  yix-xiy  f2(x-y)  =  y2x-x2y 


gi(x.y)  =  r2-x,x-y1y  g2(x.y)=  r2  -x2x-y2.v- 


(71) 


Neglecting  index  D  we  may  then  write  Eq.  (69)  in  the  form 

ft2(x-y)  f22(x-y) 

g2(x.y)  +  r2z2“  g^(x.y)  +  r2z2 

or.  after  some  reductions 

(figa  +  f2gi  )(fig2  ’  f2gi )  +  r2 (fia*i  -  \A)  =  0. 

Next  it  appears  from  Eqs.  (5)  and  (6)  that 
f,g2  1-  f2gi  =  -[Ax2  -  Ay2  +2Cxy  +  rz(2Dx  +  2Ey)j 


(72) 


(73) 


(74) 


while 


figi  -  f2gi  =  r2(2Fx  +  2Gy  +  K) 


(75) 


where 


F  =  0.5(y,  -y2):  G  =  -0.5(x,  - x2):  K  =  x,y2-x2y, 


(76) 
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Furthermore  the  last  term  of  the  equation  may  be  written 


r2(f2z2  -  f2z2)  =  -r2(Lx2  +-2Mxy  +  Ny2) 


(77) 


where 


L  =  Y2Z?  -y?Z2 


2„2 


=  Xiy,z|  -  x2y2z? 


N  =  x|z2  -  x|z2 


such  that  Eq.  (72)  may  be  written 

(Ax2  -  Ay2  +  2Cxy  +  r2(2Dx  +  2Ey))(2Fx  +  2Gy  +  K)  +  Lx2  +  Mxy  +  Ny2  =0.  j78j 

It  thus  appears  that  the  curves,  which  in  Alhazens  problem  were  conic  sections,  in  the  case 
of  diffraction  are  algebraic  curves  of  the  third  degree.  Performing  the  multiplications  and 
rearranging  the  terms  we  obtain 

a,x3+a2x2y  +  a3xy2  +a4y3  +  a5x2  +a6xy  +  a7y2 

+  r2(aax2  +a9xy  +  a10y2  +allx  +  a12y)  =  0,  (79) 


where  the  coefficients  are  given  by 


a,  =  2AF 

a4  =  -2  AG 
a7  =-AK  +  N 
a10  =  4GE 


a2 =  4CF  +  2 AG 

a5  =  AK  +  L 
a8  =  4DF 
a,,  =2DK 


a3  -  — 2 AF  4GC 

a6  =  2CK  +2M 
a9  =  4EF  +  4DG 
a,2  =  2EK. 


(80) 


where  the  various  constants  are  given  by  Eqs.  (7),  (74)  and  (77). 
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4.2  The  Alhazen  Curve 


Rearranging  the  terms  of  Eq.  (79)  we  get 

b,x3  +  b2x2y  +  b3xy2  +  b4y3  +  b5x2  +  b6xy  +  b7y2  +  b8x  +  b9y  =  0 

where 

bj  =  a[  b2  =  a2 

b4=a4  b5  =  a5  +  r2a9 

b7=a7  +  r2a10  b8  =  r2au 

From  the  y-equation 
b4y3  +  (b3x  +  b7)y2  +(b2x2  +  b6x  +  b9)y  +  b,x3  +b5x2  +b8x  =  0. 


b3  -  a3 

b6  =  a6  +  r2a9 

b9  =  r2a12. 


(81) 


(82) 


it  appears  that  there  is  a  maximum  of  three  y-values  for  each  value  of  x.  There  is  no  vertical 
asymptote  and  the  curve  intersects  the  x-axis  where 

x(b)X2  +  b5x  +  b8)  =  0. 


that  is.  at  the  origin  plus  possibly  at  two  more  points.  Similarly  the  x-equation 
b]X3  +(b2v  +  b5)x2  +(b3y2  +  b6y  +  b8)x +  (b4y3  +  b5y2  +  b9y)  =  0 


shows  that  there  are  in  genera)  three  x-values  for  each  value  of  y.  There  is  no  horizontal 
asymptote  and  the  curve  intersects  the  y-axis  where 

b4y2  +  b5y  +  b9  =0. 


that  is.  apart  from  the  origin  at  two  points  at  the  most.  Finally  the  t-equation 


b,)x2  +(b7t2  +  b6t  +  b5jx  +  b9t  +  b5  =0 


(84) 


x[(b4t3  +  b3t2  +  b2t  + 

shows  that  in  general  there  are  two  curve  points  in  each  direction  and  that  the  origin  is  not  a 
double  point.  Directions  may  exist  where  the  curve  is  not  defined  and  finally  in  general  there  are 
three  oblique  asymptotes. 

The  fact  that  there  are  three  oblique  asymptotes  suggests  that  six  diffraction  points  may  be 
possible.  We  can,  however,  eliminate  two  of  these  by  the  following  argument.  The  two  points  P, 
and  P2  must  not  be  located  on  the  same  Keller  cone,  but  rather,  one  point  must  be  on  each  of  the 
two  cones  given  by  Eq.  (70).  This  implies  that  the  points  must  not  be  on  the  same  side  of  a  plane 
perpendicular  to  the  axis  of  the  Keller  cone  at  the  point  of  diffraction,  which  means  that  their 
projections  (xlt  yj  and  (x2,  y2)  must  not  be  on  the  same  side  of  the  line  OD.  Assuming  coordinates 
(xD.  Vq)  for  the  diffraction  points  leads  to  the  condition 


(yD*i  -xDy1)(yDx2  -xDy2)<0. 


(85) 


The  condition  may  also  be  formulated  as  follows:  A  diffraction  point  can  never  be  on  the  part 
of  the  circular  edge  that  corresponds  to  the  neighboring  angles  to  the  angle  defined  by  the  origin 
and  the  radii  to  the  projections  of  P,  and  P2  on  the  plane  of  the  edge. 

Figure  1 1  shows  an  example  of  Alhazen's  curve  for  (x,.  y,.  z,)  =  (5.  15.  10),  (x2.  y2.  Z2)  =  (-5. 
-5.  3)  and  a  circle  with  r  =  10.  It  is  noted  that  there  are  six  points  of  intersection  with  the  circle. 

In  Figure  12  the  condition  Eq.  (85)  has  been  applied.  Two  of  the  intersection  points  are  now 
excluded  and  the  remaining  four  points  correspond  to  the  diffraction  points  found  by  a  numerical 
search  routine. 
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DIFFRACTION  BY  CIRCLE:  INPUT  DATA 


x,  =  5.000 

y,  =  15.000  z. 

= 

10.000 

x2  =  -5.000 

y2  =  -5.000  z2 

= 

3.000 

DIFFRACTION  POINTS  AT  FI  =  82.81 

COEFFICIENTS  OF  POLYNOMIALS 

178.06  208.86 

313.42 

a  =  -0.10000E+03 

c  =  -0.50000E+02 

d  = 

-0.50000E+01 

e  =  0.00000E+00 

f  =  0.10000E+02 

g  =  -0.50000E+01 

k  = 

0.50000E+02 

1  =  0.47500E+03 

m  =  -0.1 8250E+04 

n  = 

0.22750E+04 

a,  =  -0.20000E+04 

a2=  -0.1 0000E+04 

a3  = 

0.30000E+04 

a4=  -0.10000E+04 

a5=  -0.45250E+04 

a$=  -0.86500E+04 

a7  = 

0.72750E+04 

a8  =  -0.20000E+03 

a9  =  0.10000E+03 

a10  =  O.OOOOOE+OO 

at1  = 

-0.50000E+03 

a12=  0.00000E+00 

b,  =  -0.20000E+04 

b2=  -C.10000E+04 

b3  = 

0.30000E+04 

b4=  -0.10000E+04 

b5=  -0.24525E+05 

b6  = 

0.13500E+04 

b7  =  0.72750E+04 

b8  =  -0.50000E+05 

b9  = 

0.00000E+00 

Figure  1 1 .  Diffraction  by  a  circular  edge.  Alhazen's  curve  and  the  locus 


36 


DIFFRACTION  BY  CIRCLE:  INPUT  DATA 


x,  =  5.000 

y,  =  15.000 

Zl  = 

10.000 

x2  =  -5.000 

y2  =  -5.000 

z2  = 

3.000 

DIFFRACTION  POINTS  AT  FI  = 

82.81 

178.06  208.86 

313.42 

COEFFICIENTS  OF  POLYNOMIALS 

a  =  -0.10000E+03 

c  =  -0.50000E+02 

d  = 

-0.50000E+01 

e  = 

0.00000E+00 

f=  0.10000E+02 

g  =  -0.50000E+01 

k  = 

0.50000E+02 

1  =  0.47500E+03 

m  =  -0.18250E+04 

n  = 

0.22750E+04 

a,  =  -0.20000E+04 

a2=  — 0.10000E+04 

a3  = 

0.30000E+04 

34  = 

-0.10000E+04 

.  as  =  -0.45250E+04 

ag=  -0.86500E+04 

a7  = 

0.72750E+04 

a8  = 

-0.20000E+03 

ag=  0.10000E+03 

a10=  0.00000E+00 

an  = 

-0.50000E+03 

a12- 

0.00000E+00 

b,  =  -0.20000E+04 

b2=  -0.10000E+04 

b3  = 

0.30000E+04 

b4=  -0.10000E+04 

b5  =  -0.24525E+05 

b6  = 

0.13500E+04 

b;  =  0.72750E+04 

b8  =— 0.50000E+05 

bg  = 

0.00000E+00 

Figure  12.  Alhazen's  curve  and  locus  with  separation  condition. 
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4.3  The  Locus 

Finally  we  consider  the  locus  for  diffraction  by  a  circular  edge.  Substituting  x2  +  y2  for  r2  in 
Eq.  (79)  and  rearranging  the  terms  we  find 


csY4  +(C4X  +  C9)y3  +(c3x2  +  c8x  -+-  C,2  )y2 


+  (c2x3  +  c7x2  +C, 

,x)y  +  c,x4  +c6x3 

+  c10x2  =0 

where 

c,  =  a8 

c2  =  a9 

^3  =  a8  +  a10 

c4  =a9 

c5  =  a10 

C6  =  ai  +an 

C7  =  a2  +a12 

c8  =a3  +all 

c9  =  a4  +a12 

C10  =  a5 

C11  “  a6 

C12  =  a7  • 

(86) 


(87) 


Eq.  (86)  contains  all  terms  of  a  general  fourth  degree  polynomial  in  two  variables  except  the  first 
degree  terms  and  the  constant  term,  which  implies  that  the  curve  contains  the  origin.  From  the 
y-equation 

c5y4  +  (c4x  +  c9)y3  +(c3x2  +c8x  +  cl2)y2 

+  (c2x3  +  c7x2+cux)y  +  c1x4+c6x3  +  c10x2  =0  (88) 

it  appears  that  there  is  a  maximum  of  four  y-values  for  each  value  of  x,  that  no  vertical  asymptote 
is  present  and  that  the  curve  intersects  the  x-axis  where 

x2(c!x2  +c6x  +  c10)  =  0. 

that  is,  at  the  origin  plus  possibly  two  more  points.  Similarly  the  x-equation 
c,x4  +(c2y  +  c6)x3  +(c3y2  +c7y +  c10)x2 

+  (c4y3  +  c8y2  +  C,  ,y)x  +  c5y4  +  c9y3  +  c12y2  =  0  (89) 

shows  that  there  are  in  general  four  x-values  for  each  value  of  y.  There  is  no  horizontal  asymptote 
and  the  curve  intersects  the  y-axis  where 
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c5y2+c9y+ci2  =o. 


(90) 


that  is.  apart  from  the  origin  at  two  points  at  the  most.  Finally  the  t-equation 

x2[(c5t4  +c4t3  +  c3t2  +c2t  +  Cj)x2  +(c9t3  +c8t2  +c7t  +  c6)x  +  e12t2  +cnt  +  c10  j  =  0  j90) 


Indicates  that  there  are  in  general  only  two  curve  points  in  each  direction.  Furthermore  the  origin 
is  a  double  point  and  a  possibility  of  four  oblique  asymptotes  exists.  Apart  from  these  very  general 
statements  it  is  difficult  to  be  specific  about  properties  of  the  curve,  except  for  the  remark  that 
the  presence  of  four  asymptotes  seems  incompatible  with  the  fact  that  there  are  only  two  curve 
points  in  each  direction. 

To  illustrate  the  locus  in  a  specific  case  it  is  shown  in  Figure  1 1  together  with  Alhazen's 
curve.  The  locus  intersects  the  circle  at  the  same  six  points  as  does  Alhazen's  curve  as  should  be 
expected.  Since  the  points  on  the  locus  are  subject  to  the  same  separation  condition  as  the  points 
of  Alhazen's  curve,  we  may  apply  the  condition  Eq.  (85)  to  the  locus  also.  The  result  is  shown  in 
Figure  12  where  only  four  diffraction  points  remain.  It  is  noted  that  the  separation  condition 
removes  an  entire  branch  of  the  locus  and  leaves  a  continuous  curve  defined  solely  by  the 
coordinates  of  PI  and  P2.  In  this  case  the  curve  indicates  the  radius,  which  separates  the  circles 
with  two  and  four  diffraction  points.  It  is  an  interesting  fact,  that  the  locus  may  split  into  two 
separate  curves,  one  of  these  being  a  closed  curve.  This  is  illustrated  in  Figure  13.  where  r  =  10, 
(x,.  y,.  z,)  =  (5.  15,  10)  and  (x2.  y2.  24)  =  (-5.  -5.  4.75).  The  closed  curve  defines  two  circles  with 
radii  rj  and  r2,  where  r,  <  r2.  Thus  for  the  given  location  of  P,  and  P2  there  will  be  two  points  of 
diffraction  if  the  radius  of  the  circle  r  >  r2  or  if  r  >  rj.  If  r*  <  r  <  r2  there  will  be  four  points  of 
diffraction. 
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DIFFRACTION  BY  CIRCLE:  INPUT  DATA 


x,  =  5.000 

y,  =  15.000 

Zl  = 

10.000 

x2  =  -5.000 

y2  =  — 5.000 

z2  = 

4.750 

DIFFRACTION  POINTS  AT  FI  = 

COEFFICIENTS  OF  POLYNOMIALS 

82.81 

311.28 

a  =  -0.10000E+03 

c  =  -0.50000E+02 

d  = 

-0.50000E+01 

e  =  O.OOOOOE+OO 

f=  0.10000E+02 

g  =  -0.50000E+01 

K  = 

0.50000E+02 

1  =  -0.25766E+04 

m  =  -0.80700E+03 

n  = 

0.19350E+04 

a,  =  -0.20000E+04 

a2=  — 0.10000E+04 

a3  = 

0.30000E+04 

a4=  -0.10000E+04 

a5  =  -0.75766E+04 

ae=  -0.661 40E+04 

a7  = 

0.69350E+04 

as  =  -0.20000E+03 

a9  =  0.10000E+03 

a,0  =  0.00000E+00 

3i,  = 

-0.50000E+03 

a12-  0.00000E+00 

b,  =  -0.20000E+04 

b2=  -0.10000E+04 

b3  = 

0.30000E+04 

b„=  -0.10000E+04 

b5  =  -0.27577E+05 

b6  = 

0.33860E+04 

b7=  0.69350E+04 

b8  =  -0.50000E+05 

b9  = 

O.OOOOOE+OO 

Figure  13.  Diffraction  by  a  circular  edge.  Separation  condition  applied.  The  curve  shown  is  the 
locus. 
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5.  CONCLUSION 

In  the  present  report  the  results  of  a  preliminary  investigation  of  some  ray  tracing  problems 
related  to  circles  and  ellipses  are  described.  The  work  is  based  on  an  alternative  formulation  of 
the  law  of  reflection,  that  leads  directly  to  an  analytic  solution  to  the  problem  of  finding  the  points 
of  reflection  on  the  concave  side  of  a  spherical  mirror.  It  is  shown  that  the  solution,  referred  to  as 
Alhazen’s  solution,  also  applies  to  reflection  by  the  convex  side  of  a  spherical  mirror. 

Furthermore  a  new  concept  in  ray  tracing,  the  locus,  is  introduced  and  applied  to  the  circle 
problem.  This  leads  to  an  alternative  rule  for  the  a  priori  determination  of  the  number  of 
reflection  points  and  suggests  a  new  numerical  method  for  the  determination  of  reflection  points. 

The  techniques  described  above  are  applied  to  the  ray  tracing  problems  associated  with 
reflection  by  an  ellipse  and  diffraction  by  a  circular  edge.  The  solution  of  tnese  problems  is  closely 
related  to  the  geometry  of  plane  algebraic  curves  given  by  polynomials  in  two  variables.  Table  1 
shows  the  occurrence  of  the  terms  in  the  polynomials  for  the  three  cases  considered,  namely 
reflection  by  a  circle  and  by  an  ellipse,  and  diffraction  by  a  circular  edge.  It  is  seen  that  the 
governing  equations  become  more  complicated  for  the  ellipse  than  for  the  circle,  and  that  the 
diffraction  problem  is  more  involved  than  the  reflection  problems.  This  illustrates  the  fact  that  in 
the  present  cases,  reflection  may,  from  a  ray  tracing  point  of  view,  be  considered  a  special  case  of 
diffraction. 

For  reflection  by  an  ellipse,  the  locus  method  is  modified  and  it  is  shown  that  no  more  than 
four  reflection  points  are  possible.  In  the  diffraction  case  formally  six  solutions  result,  but  two  of 
these  may  be  eliminated,  since  they  correspond  to  non-physical  configurations.  Application  of 
this  observation  may  on  the  average  halve  the  computer  time  spent  on  ray  tracing  in  this  case. 

Further  research  is  needed  to  assess  the  significance  of  the  results  found  so  far.  Topics  for 
further  work  include: 

1.  Simplifications  resulting  from  the  field  point  being  at  infinity. 

2.  Ray  tracing  for  diffraction  by  an  ellipse. 

3.  Multiple  reflection  and  diffraction. 

4.  Application  of  an  efficient  numerical  technique  to  the  locus  method. 

5.  Generalization  to  three  dimensional  problems. 
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